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I. INTRODUCTION 


Without comparison, the Standard Model of particle physics (SM) is the most success¬ 
ful physics model to date, accurately predicting an enormous number of observables with 
high precision from only a handful of htted parameters. The success of the SM may have 
culminated in 2012 when the ATLAS and CMS experiments announced the discovery of 
the Higgs boson [l, 2|, predicted by the SM as a direct result of the electroweak symmetry 
breaking which was introduced to provide masses into the theory. Still, there are a number 
of observations which may not be explained within the SM itself. Most notable among these 
are the existence of dark matter, the exclusion of gravity and the observation of neutrino 
oscillations. In addition, there are conceptual theoretical problems with the SM, such as 
the hierarchy problem, indicating that the SM may only be a low-energy approximation of 
a more general theory. As such, the SM should be viewed as an effective theory and a priori 
higher-dimensional operators, suppressed by powers of a new mass scale A should be added 
to the SM Lagrangian. At lower energies, the additional effective operators will generally 
produce very small corrections due to this suppression. This concept is further supported 
by the fact that the only gauge invariant operator allowed at dimension hve, and therefore 
suppressed only by one power of A, is the so-called Weinberg operator js], which results in 
a Majorana mass term for the left-handed neutrinos of the SM. It is therefore not unrea¬ 
sonable to imagine that the effect of neutrino masses would be among the hrst observations 
of physics beyond the SM, which indeed is the case due to neutrino oscillations requiring 
neutrino mass-squared differences to be non-zero. 

Neutrino flavour conversion, although at that time not conhrmed as such, was hrst ob¬ 
served in solar neutrino experiments where a discrepancy between the observed hux and the 
hux predicted by solar models was found [^. Since the hrst robust evidence of neutrino 
oscillations by the Super-Kamiokande experiment’s observation of atmospheric neutrinos in 
1998 j^, they have been extensively studied experimentally in a variety of atmospheric, so¬ 
lar, reactor, and accelerator experiments, which have helped to constrain the neutrino mass 
and mixing parameters to very high precision (see Refs, bijs] for recent global hts). The 
remaining questions in neutrino oscillation physics today are the neutrino mass ordering, 
the existence or non-existence of CP violation in the lepton sector and the octant of the 
leptonic mixing angle 6 ^ 23 . Answering these three questions is the main aim of the next gen- 
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eration of neutrino oscillation experiments, such as the European Spallation Source Neutrino 
Super-Beam (ESSz/SB) experiment {qI, which is a proposed accelerator neutrino experiment 
based on the European Spallation Source (ESS) currently under construction in Lund, Swe¬ 
den. The sensitivity of ESSi^SB to the CP-violating phase 6 was studied in Ref. 1^, while 
the sensitivity to other standard oscillation parameters was discussed in Ref. and the 


sensitivity to light sterile neutrinos in Ref. 11|. 


While the Weinberg operator provides the neutrino masses necessary for neutrino os¬ 
cillations to occur and neutrino oscillations have been hrmly established as the leading 
mechanism behind neutrino flavour conversion, higher-order operators may give rise to sub¬ 
leading contributions to the neutrino conversion probabilities and their observation would 
allow us to gain additional insight into the high-energy completion of the SM and the gener¬ 
ation of neutrino masses. In addition, it may be necessary to consider the robustness of the 
usual neutrino oscillation parameters when higher-order operators are also considered. One 
of the more common types of operators to be investigated in this respect is non-standard 
neutrino interactions (NSIs), which are effective four-fermion operators involving at least 
one neutrino held. For recent reviews on NSIs, see Refs. [^, I^ . 


In this work, we will consider the possible impact of NSIs at the ESSi^SB experiment. 
We will study both the inhuence of NSIs on the determination of the standard neutrino 
oscillation parameters and the bounds which ESSz/SB could place on the NSI parameters. 
In particular, we will focus on correlations in the determination of the leptonic CP violation 
and the NSI parameters, which is of large importance for ESSz/SB as the discovery of leptonic 
CP violation is the major scientihc target of this experiment. 

The rest of this work is organized as follows. In Sec. m we will briehy review non¬ 
standard neutrino interactions and present the current upper bounds on the source and 
detector NSI parameters. Next, in Sec. IIIIl the setup of the proposed ESSz/SB experiment 
will be discussed. Then, in Sec. lYl we will investigate the phenomenology of source and 
detector NSIs at probability and event-rate levels. In Sec. IVj the main results of our full 
computation on source and detector NSIs at ESSz/SB will be presented. Finally, in Sec. m 
we will summarize and draw our conclusions. 
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II. NON-STANDARD NEUTRINO INTERACTIONS 


When considering NSIs, we will be confronted with effective four-fermion operators of 
the type 

= {fll^PL,Rf2)if3l,PL,Rh) + h.C. , (1) 

where fi {i = 1,2, 3,4) are SM fermion helds and Pl,r are left- and right-handed projec¬ 
tions. These operators are of dimension six and they will therefore appear together with an 
effective coupling constant of dimension minus two in the effective Lagrangian. Since we are 
interested in the NSIs of neutrinos, we require that at least one of the fermion helds in the 
operators is a neutrino held, which implies that the corresponding projection operator must 
be Pl- Furthermore, in order to keep the electromagnetic and strong interactions unbroken, 
we require that all operators are scalars under transformations of the corresponding gauge 
groups. Due to the weak interaction being broken, we do not impose any constraints on 
the transformation of the operators under SU(2 )l. It should be mentioned that imposing 
SU(2 )l gauge symmetry on the dimension-six operators would lead to havour constraints on 
these operators 1^, ^ , leaving only a few possible operators without signihcant constraints 
due to the non-observation of ehective four-charged-fermion processes such as /i —?• 3e. 
The dimension-six operators which break 811(2)^ may generally be induced from higher¬ 
dimensional operators such as (00^)0, where 0 is the Higgs held, which are invariant under 
SU(2 )l, but generate SU(2)L-breaking terms once the Higgs held takes on a vacuum expec¬ 
tation value V. Depending on the dimension at which the NSIs are generated above the 
electroweak scale, we may expect the NSI coefficients to scale as where n is the 

dimension and A is the energy scale at which the NSIs are generated . 


The diherent possible neutrino NSIs are generally c 
four-fermion operators. The neutral-current NSIs la. 


ivided into two categories of ehective 

17 | 


(9 


f{L,R) _ 
aj3 


= {jya'y^PLJ^fi){f'ytiPL,Rf) + h.c. , 


( 2 ) 


where / is a charged fermion held, ahects the neutrino havour propagation in matter for / = 
u, d, e by providin g a n ehective potential analogous to the Mikheyev-Smirnov-Wolfenstein 


(MSW) potential 16 


18 


19| . For the neutral-current NSIs to be of importance, relatively 


large matter potentials and/or high neutrino energies are required. As this is not the case 
for the ESSz/SB experiment, we will not focus on such NSIs in this work. On the other hand. 
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the charged-current NS Is 201 




( 3 ) 


where / and f are different fermion helds such that the operator is invariant under U(l)gf^ 
and SU(3)^, will instead affect the production and detection processes of neutrinos and this 
effect will not depend on the neutrino energy or the presence of matter along the neutrino 
propagation. 

In the remainder of this work, we will focus on the charged-current NSI Lagrangian 

£nsi = -2V2Gf Y.^tp{~Pl^Pda){d^,Pxu) + h.c. , (4) 

X&{L,R} a, 13 

which includes the operators that will appear in neutrino production by pion decays vr —)■ £q,z/ 
and charged-current neutrino detection processes. Here, we have normalised the strength of 
the NSIs to that of the weak interaction by the introduction of the Fermi coupling constant 
Gp. The NSI parameters Sais are therefore dimensionless numbers expected to be of the 
order (u/A)"’“^. With the introduction of charged-current NSIs, the production amplitude 
of the neutrino mass eigenstate |z/j) in the vr’*', which in the SM is proportional to 17*^, where 
U is the leptonic mixing matrix, is now instead proportional to where 

the NSI parameters relevant for the source process are 


^q(3 


^^af3 ■ ( 5 ) 

Unlike the source process, the detection process does not necessarily involve a pseudoscalar 
current in the quark sector. We instead dehne the NSI parameters relevant for the detection 
process as 


eU = iV > ( 6 ) 

where P represents the quark current in the detection process. Due to the nature of the 
inverse beta decay involved in the detection process, this dehnition oversimplihes the neu¬ 
trino oscillation probabilities that we will discuss in Sec. lIVl^l However, we will use this as 


^ In fact, the neutrino oscillation probabilities should be computed along the lines 

Pap ^ ^ [Pap{e^,e^) -b4.5Pa/3(e^,e^)] , 

where Pap{s^, e'^) is the probability for a given source/detector NSI. Note that the largest prefactor comes 
from the contribution with the source and detector effects both dependent on the axial quark current. 
This would therefore indicate a relation between the source and detector NSIs. 
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a simplified model for how NSIs may affect ESSz/SB. The complex conjugate and change of 
indices has been introduced to adhere to the usual convention in the held when considering 
detector NSI effects. The production rates of charged leptons at the detector in any neu¬ 
trino oscillation experiment will be affected by this change in the production and detection 
amplitudes and we may ask the question whether or not the presence of such NSIs could be 
measured or have a negative impact on the experimental precision to the standard oscillation 
parameters. The experimental bounds (at the 90% C.L.) on the NSI parameters relevant 
for the ESSz/SB experiment from non-oscillation experiments are given by 211 


< 0.026 , KJ < 0.078 , \eU < 0.013 , 


-.d 


s 

flT\ 


\ei\ < 0.041 , < 0.025 , |4,| < 0.041 , 

< 0.026 , < 0.078 , \etJ < 0.013 . 


(7) 


Although these bounds are quite stringent, it should be kept in mind that the next gener¬ 
ation of neutrino experiments is aiming for highly sensitive measurements of the neutrino 
oscillation parameters. As such, even sub-leading effects may be of interest and it is worth 
the effort to examine the possible impact of these effects. It is also worth noting that new 
oscillation experiments, such as those performed with nuclear reactors, may be sensitive 


to some of these NSI parameters as well 


22l |. However, the current bounds from these 


experiments are somewhat weaker than the bounds quoted above 231 


III. THE ESSz/SB EXPERIMENTAL SETUP 


In this section, we describe the experimental setup for the proposed ESSz/SB experiment. 
We have used the standard flux (with 2 GeV protons) and cross-sections from the ESSz/SB 
collaboration {q]. The source provides a neutrino beam for two years and an antineutrino 
beam for eight years. We have assumed that a 500 kiloton water Cherenkov detector is placed 
at a distance of 540 km from the source, which corresponds to the location of the mine in 
Garpenberg, Sweden. The detector specihcations have been taken from the performance 


study of the MEMPHYS detector 


24l |. The energy range of interest is up to 2 GeV, which 


is divided into 20 energy bins. We have used 9% (18%) systematic errors on the signal 
(background) events. Unless specihed otherwise, we have also assumed the existence of a 
near detector with mass 1 kiloton, 1 km from the source and the same flux as at the detector 
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at 540 km scaled by the distance-squared. As a crude approximation, we assume the same 
characteristics for both these detectors. 

To this end, we have written our own probability engine to calculate the neutrino os¬ 
cillation probability in the presence of source and detector NSIs. This probability engine 


interfaces with GLoBES 


25 


26| for calculating the neutrino event rates at ESSi^SB. The 


large parameter space is handled with the help of the GLoBES plugin MonteGUBES [27|. 


IV. NEUTRINO OSCILLATIONS WITH NSIS 


Standard three-flavour neutrino oscillations depend on six fundamental parameters - two 
mass-squared differences, and Amh, three mixing angles, 6 ^ 12 , 6^13 and 623 and one 

GP-violating phase 6 . In addition, if the neutrinos are propagating through matter, the 
charged-current interactions of the neutrinos with electrons modify the oscillations. This 
effect can be incorporated into the probability formalism by using the MSW potential term 


A = 2\/2Gpn^E ly, llSl, ll9|, where Ug is the number density of electrons in the matter and 


E is the neutrino energy. For an experiment like ESSz/SB with a short baseline length as 
well as low neutrino energy, we can ignore the matter effects for the sake of this discussion. 
(The numerical results presented in this work do not make any such assumption.) 

Non-standard neutrino interactions can affect the production and detection of neutrinos 
at the source and detector, respectively. In the SM, interactions of charged leptons with 
neutrinos are strictly flavour-diagonal. However, charged-current NSIs can introduce a non¬ 
zero overlap between charged leptons and neutrinos of different flavours. Thus, a neutrino 
produced at a source in association with a charged lepton ia is not simply but is given 
by 


20 


28 


30| 


\0 = Wa) + 


'a 7 1 7 




7=e,/i,T 


Similarly, a neutrino that produces a charged lepton in a detector is 


(z/^i — (z/^i + 


^7/3 W71 


( 8 ) 


(9) 


7=e,/L,T 


The matrices e* and are in general complex, giving 36 new parameters. These are 9 
amplitudes and 9 phases of each NSI parameter in the source and detector NSI matrices. 

Not all of the 36 NSI parameters are relevant for the experiment under consideration. 
Since we are only interested in the oscillation channels —)■ and z/^ —)■ (and their GP 
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conjugates), the relevant parameters are and where 7 G {e,p,r}. Thus, the 

parameter space is reduced to 9 complex or 18 real parameters, in addition to the standard 
ones. In this work, we treat all of them as independent parameters, which is the most general 


case. 


Deriving an analytical formula for the neutrino oscillation probabilities is difficult even in 
the standard three-flavour scenario. Typically, expressions for the probabilities are given as 
perturbative expansions in small parameters such as or sin 6*13 3ll-l33|. For the 

discussion in this section, we refer to the analytical formulae for vacuum oscillation probabil¬ 
ities derived in Ref. 3^, which include source and detector NSIs. These formulae are valid 
up to second order in Am^i/Amg]^ and sin 6 ^ 13 , and up to hrst order in the NSI parameters. It 
is easy to observe that linearizing the expressions in the NSI parameters and ignoring cubic 
and higher order terms overall, leaves only a few NSI parameters in the expressions. For 
instance, in the case of the vacuum probability only the NSI parameters and £^g 

are present up to linear order. While these approximate analytical formulae provide useful 
insights into the physics of NSIs in neutrino oscillations, we stress that all simulation results 
presented in this work make use of numerically computed neutrino oscillation probabilities 
without approximations. 

In Fig.lll we have plotted the variation of the neutrino oscillation probability P^e with the 
amplitude of each of the relevant NSI parameters. The range of values chosen for the NSI 
parameters is the 90% C.L. bounds on them as listed in Eq. ([7]). Each of the probabilities 
shown are calculated numerically, using <5 = 0, ^23 = 45° and normal neutrino mass ordering; 
and all other NSI parameters, including phases, set to zero. Out of the three NSI parameters 
present up to linear order, the variation due to e^g is the strongest, while that due to 
is the weakest. This pattern follows from the allowed range given by the current bounds. 
Out of the remaining three, has the greatest effect, which is again because it is not very 
tightly constrained by current data. 

Figure [T] is plotted for a hxed value of 5 = 0. For ESSz/SB, it is interesting to explore the 
interplay between 5 and the NSI parameters. To this end, we show in Fig. [2] bi-probability 
plots for ESSi/SB. This hgure is presented for a hxed energy of 400 MeV, which corresponds 
to the second oscillation maximum for the ESSi^SB baseline. This is also the energy around 
which the unoscillated event rate is maximal. As 5 varies over its full range, the neutrino 
and antineutrino probabilities trace out an ellipse as shown. In the standard case, we obtain 











0.2 0.4 0.6 0.8 1 
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FIG. 1. Neutrino oscillation probability as a function of the neutrino energy E and its variation 
with each of the relevant NSI parameters. The values of the NSI parameters are chosen within their 90% 
C.L. bounds, i.e. assuming their phases to be either 0 or tt. The variation is shown for both neutrinos 
and antineutrinos. The values of the fundamental neutrino parameters are set to = 7-6 x 10“®eV^, 

Am§i = 2.45 X eV^ sin^ 9^2 = 0.304, 6»23 = 45°, sin^ 26»i3 = 0.09 and 5 = 0. 


the central (blue) ellipse. In each of the panels of this hgure, one NSI parameter is varied 
within its 90% C.L. bound, which gives the spread in the ellipse. 

In order to explain the features observed in Fig. O we dehne the variation of the neutrino 
oscillation probability as 


- P'r(el, = 0 ) , 


( 10 ) 


where a,(3 E {e,/t, r} and x G Using the perturbative analytical expression for P™ 
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from Ref. [3^, we obtain for the cases of and 


- -4|£;,g| sin 013 sin 023 sin (A+ 5) sin A , 
^^;r(<e) - -M4e\ sin ^13 cos 2023 sm023 COS(5sin^ A 


—, 


( 11 ) 


Ap;r(4.) - 



( 12 ) 


( 13 ) 


where A = AmliL/{4E). In deriving each of Eqs. ffTT]) - ffT^ . we have set all other NSI 
parameters to zero. Note that for the cases of and efg, there are no linear-order 

terms in the corresponding formulae, and the dependence on the NSI parameters only appear 
at second order and above. First, we observe (as in Fig. [T]) that the variation of P^g is the 


largest for £^g (due to linear variation and weakest upper bound) and the smallest for 


(due to higher-order variation and strongest upper bound). For e®g, efg and £((g, the 
variations are intermediate, depending on a non-trivial combination between the value of 
the upper bound on the considered NSI parameter and if this NSI parameter appears at 
linear order or not in the variation. Second, we can explain the structure of the band for 
each panel. We illustrate this for the case of £^g. For the baseline and energy considered, A 
evaluates to around —120°, close to the second oscillation maximum. It is then easy to see 
that the maximum ‘width’ of the band occurs when A -|- 5 = ±90°, i.e. when S is around 30° 
or —150°. Likewise, for A ± 5 = 0,180°, the probability becomes independent of £®g, and 
the band ‘pinches off’. This occurs when 6 is around 120° or —60°. For antineutrinos, the 
sign of 6 is changed, and one can use similar arguments to hnd the broadest and narrowest 
points along the P^g axis as well. 

Finally, Fig. [3] shows the neutrino and antineutrino event rates for FSSi^SB in a bi-rate 
plot, using the same parameter values as for Fig. |2l The event rates plotted are the total 
rates across all energy bins. Statistical error bars have been included for four representative 
values of 6 . In addition to the variation of the probabilities in Fig. |2l this hgure gives a hrst 
indication of the impact of NSIs versus the possible experimental resolution of the FSSi^SB. 
Where the experimental error bars on the total event rates are smaller than the possible 
variation of the NSI parameters, the ESSi^SB will generally be sensitive to NSIs smaller 
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FIG. 2. Bi-probability ellipse for ESS:/SB and its variation with the relevant NSI parameters. The values 
of the fundamental neutrino parameters and the NSI parameters are the same as in Fig. [1] 


than the current bounds. However, note that the converse is not necessarily true as the 
experimental results do not only include the total event rates, but also spectral information, 
which may also be used to constrain the NSIs. In particular, this will be apparent for our 
results on which does not change the event rates signihcantly. 


V. RESULTS ON NSIS AT ESSnSB 


The main goal of the proposed ESSz/SB experiment is to measure the CP-violating phase 
6 with high precision. In this section, we examine both the impact of the NSI parameters 
on this 6 measurement and the ability of ESSi^SB to measure the NSI parameters. 

The central values of the neutrino parameters Am^i, |Am|;^|, 612 and 615 are taken close to 


their current best-£t values 


a-jsl . We have also imposed Gaussian priors on these parameters 


with a width obtained from these global hts. The values of 623 and 6 used are different in 
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FIG. 3. Bi-rates ellipse for ESS:/SB and its variation with the relevant NSI parameters. The values of the 
fundamental neutrino parameters and the NSI parameters are the same as in Fig. [TJ 


each case, and are specified in the text. In addition, we have assumed a 5% prior on the 
true value of sin^ 26 * 23 . The NSI parameters are of the form where a,/9 G {e,/i,r} and 
X G {s,d}, since the source and detector NSI parameters can be different in general. Thus, 
we have 18 complex NSI parameters, or 36 real NSI parameters, in addition to the standard 
ones. We have run our simulations for both normal (NO) and inverted (10) neutrino mass 
ordering. We find that there is very little qualitative difference between the results in these 
two cases. Therefore, in what follows, we show only the NO results. 


A. Effect on precision measnrement at ESSi^SB 

In this subsection, we discuss the interplay between NSI parameters and the 6 precision 
of ESSz/SB. The results are shown in the form of precision contours in the O 23 — S plane. This 
is performed for three representative values of 623 G {42°, 45°, 48°}; and four representative 
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values of 5 G {-90°, 0°, 90°, 180°}. 

First, we explore the effect of marginalizing over the source and detector NSI parameters 
on precision measurements at ESSi^SB, in the special case the true NSI parameters are zero. 
In other words, we take all the NSI parameters to be zero when generating the mock data, 
but allow them all to vary in the £t. Thus, these plots show the robustness of the ESSz/SB 
measurements against a scan for NSIs. The results are shown in Fig. 01 The solid curves 
show the 68%, 90% and 95% C.L. contours for the allowed region in the parameter space. 
The dashed contours are for the standard case where there are no NSI parameters in the 
data or the £t. 

We observe that the search for NSIs does not affect the 623 precision of ESSz/SB much. 
The precision in 6 is worsened to at most twice its standard precision, in the worst case. 
For most cases, the precision is seen to be quite robust, even in spite of a severely enlarged 
parameter space. This is true, irrespective of the true value of 623 or 6 . 

Second, in Fig. O we investigate the same effect as in Fig. 01 but with a non-zero value 
for the NSI parameters in the mock data. These ‘true’ values of the NSI parameters have 
been taken to be half of the bounds given in Eq. ([7]) for the amplitudes. The true values 
of the non-standard phases are taken to be zero. In the £t, as before, all the standard as 
well as the NSI parameters are marginalized over. Thus, these plots show the robustness of 
measurements at ESSz/SB against a scan for NSI parameters, but in the presence of NSIs. 

As in the previous case, in the presence of NSIs, we observe that the 623 measurement is 
not affected much, while the precision in 6 worsens. Here, the worsening depends signihcantly 
on the value of 6 in nature. When 5 = 0, the worsening of precision is least, whereas for 
6 = 180°, the precision is worst. This is seen uniformly across the range of 623 values 
considered. The reason for this is as follows. A measure of the precision of 6 is dP™°/d5. In 
order to hnd the value of 6 at which this precision is minimal, we set the derivative of this 
quantity, i.e. d^P™°/d5^ to zero. Since the dependence of the probability on 5 is harmonic, 
the second derivative is proportional to the probability itself. As seen from the panels in 
Fig. O the smallest probability for both neutrinos and antineutrinos is around 180°. This is 
why the precision of 6 is worst at 180° in the presence of NSIs. 

Third, we study how the precision measurement at ESSz/SB would be affected if NSIs are 
present in nature, but are not accounted for in the scan of the parameter space. For this, 
we have taken non-zero values of the NSI parameters in the mock data (the same non-zero 
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values as in the previous case of Fig. |5]), but their values have been kept hxed at zero in 
the £t. The results are displayed in Fig. [ 6 l Here, the solid curves represent the 68 %, 90% 
and 95% C.L. contours, when NSIs are present in the data, but not in the £t. The dashed 
curves are the corresponding contours for the case where the NSIs are marginalized over in 
the £t. Thus, the difference between the solid and dashed contours indicates the effect of our 
ignorance of the existence of NSIs. Our ignorance leads us to an over-optimistic precision 
in 6 , as expected. The effect is more pronounced for the true value of 5 = 180°. As before, 
the 6*23 precision is not affected. 


B. Constraining NSI parameters at ESS 


Having investigated the effect of NSIs on precision measurements at ESSi^SB, we explore 
the ability of this experiment to measure the NSI parameters themselves. As we have seen 
before, the effect of the NSI parameters on the probability is quite mild. Therefore, we do 
not expect to obtain very strong constraints on these parameters. 

Figure [7] shows the limits which ESSi^SB can set on the amplitudes of the NSI parameters 
for NO, 6*23 = 45° and 5 = 0. Consider the top-left panel, corresponding to the parameter 
In generating this plot, we have set the true values of the NSI parameters to be zero. We 
show the a function of the test value of |e®e|, when all the other neutrino parameters, 
including the NSI ones are marginalized over. Horizontal lines have been drawn in the plots, 
corresponding to 68 %, 90% and 95% C.L., assuming a distribution. One can read off 
the limits that ESSz/SB can impose on these parameters from this plot. Similarly, the other 
panels show the limits for the other relevant parameters. 


The 90% C.L. limits on the NSI parameters using data from ESSi^SB are summarized 
in Table E The hrst column gives the limits when all the other NSI parameters are kept 
free in the £t, which can be simply read off from Fig. [71 These limits should be interpreted 
as being realistic, since they are derived without making any assumptions on the values of 
the other NSI parameters. We have also computed the limits when the NSI parameters are 
only considered one at a time, i.e. all other NSI parameters are hxed to zero. These limits, 
which are given in the second column, are more optimistic. The realistic limits on |e^g| and 
£^g| are comparable to the ones in Ref. 2l|], which are listed in the third column for ease 
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risonlfl 


of comparisono This is because these parameters have the maximum effect on P^e, as seen 
from the analytical expressions. In the optimistic case, the limits on |£^g| and |£„g| improve 


"fie \ 


by a factor of two compared to the existing bounds. For all the other NSI parameters (except 
|e®g| and the realistic and optimistic limits basically coincide and are less stringent 


than the limits in Ref. 


2l|. 


Parameter 

Limits with all other 

NSI parameters free 

Limits with all other 

NSI parameters zero 

Limits from Ref. 

\s^ 1 

\^fie\ 

0.025 

0.014 

0.026 

le* 1 

0.27 

0.27 

0.078 

1 

l‘-/iT I 

0.040 

0.040 

0.013 

\e^ 1 

Peel 

0.15 

0.15 

0.041 

1 ^^ 1 

Pe/il 

0.087 

0.082 

0.026 

\e^ I 

P/iel 

0.025 

0.014 

0.025 


0.28 

0.27 

0.078 

1 

Pre 1 

0.11 

0.12 

0.041 

Pr/il 

0.040 

0.033 

0.013 


TABLE I. 90% C.L. sensitivities of ESSi^SB to the NSI parameters. 


We have also checked whether ESSi^SB can measure the values of the NSI parameters 
with any reasonable precision. The procedure for this is the same as for Fig. [71 except that 
the true values of the NSI parameters are non-zero. We have chosen these true values to be 
half of the 90% C.L. bounds given in Ref. 2l|. The resulting from these computations is 
shown in Figs. |8]and|9]for the amplitudes and the phases of the NSI parameters, respectively. 
We hnd that ESSi^SB is not capable of distinguishing the chosen non-zero values of the 
parameters from zero, even at 68% C.L., nor is it able to significantly constrain any of the 
NSI phases. 


Note that the limit on a given NSI parameter in Ref. 2l| has been computed considering only that 


parameter and assuming all other NSI parameters to be zero, which corresponds to the optimistic case. 
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C. Role of near detector and systematics 


Throughout this study we have used systematic errors of 9% in the signal and 18% in 
the background events 2^. These are typical values for a superbeam experiment with a 
megaton-scale water Cherenkov detector. In order to study the role of systematic errors on 
our results, we have also simulated our experiment with a smaller systematic error of 5% 
in both signal and background. This is of course a very optimistic value. We have found 
that the limits on NSIs from ESSi^SB do not change appreciably with this drastic reduction 
of systematic errors. This is because in spite of having a large detector and intense source, 
ESSz/SB is still statistics-dominated due to the lower event rate at the second oscillation 
maximum. 

Finally, we examine the role played by the near detector in the sensitivity of ESSz/SB. 
As described before, we have used a crude simulation of the near detector throughout this 
work. Here, we compare the results of our simulation with and without the near detector. 
We show in Fig. [10] a recomputed version of Fig. [S] both with and without the near detector. 
The solid contours are the same as before, but the dashed contours show the same allowed 
regions, if only the far detector is used. We observe that in the absence of a near detector, the 
5-sensitivity of ESSz/SB is worsened. The limits on NSI parameters are also worse without 
the near detector. 


VI. SUMMARY AND CONCLUSIONS 

In this work, we have investigated the effects of source and detector NSIs at the proposed 
neutrino oscillation experiment ESSz/SB, with a baseline of 540 km - the source being 
the ESS in Lund, Sweden and a MEMPHYS-like detector in Garpenberg, Sweden. The 
ESSz/SB experiment is designed to determine the leptonic CP-violating phase 5 at the second 
oscillation maximum. However, it may also be able to probe source and detector NSIs. Due 
to the short baseline length and low neutrino energy of this experiment, matter NSIs will 
not be of importance, and are therefore not considered in this work. 

First, we have studied the three-flavour neutrino oscillation probabilities with source and 
detector NSIs, which depend on six relevant NSI parameters - efg, £^g and £^g. 

We used perturbative analytical expressions for the z/^ —)■ z/g channel that is the important 
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channel for ESSz/SB in which e^g and e^g are the dominating NSI parameters in order to 
observe the impact of these parameters. We have found that, for the range of values allowed 
by the current data, the NSI parameter e^g affect this probability the most, whereas the NSI 
parameter e^g the least. All other four NSI parameters have intermediate influence on the 
probability. 

Second, we have explored the effect of marginalizing over the NSI parameters on precision 
measurements at ESSz/SB using two cases: (i) The true values of the NSI parameters are 
set to zero and (ii) the true values are set to half of the current 90% C.L. bounds. In both 
cases, the precision of measuring S is reduced by at most a factor of two. In addition, a 
measurement of the leptonic mixing angle 623 is not affected by NSIs. If we do not take 
the effect of NSIs into account when determining the value of S, we obtain over-optimistic 
results. The effect is most pronounced for a true value of h = 180°. Note that the impact of 
NSIs on the results are qualitatively same for both NO and 10. 

Third, we have determined the possibility of ESSz/SB to measure the values of the NSI 
parameters. In a realistic case with all NSI parameters free, we have found limits on e^g 
and e^g at 90% C.L. that are similar to the existing limits in the literature, whereas in a 
optimistic case with only one NSI parameter free and the rest set to zero, the limits on 
£^g and £^g are improved by a factor of two. Furthermore, we have set the true values of 
the NSI parameters to half of their existing bound, and found that ESSz/SB is not able to 
differentiate the set values from zero at 68% C.L. or impose any signihcant constraints on 
the phases of the NSI parameters. 

Finally, we have examined the influence of the presence of a near detector at the ESSz/SB 
experimental setup. Indeed, we show that without a near detector the results would be more 
pessimistic concerning both the sensitivity of S and the limits on the NSI parameters. Note 
that the results are not changed signihcantly by reducing the systematic errors. 

In conclusion, using ESSz/SB with a near detector, the presence of NSIs will at most 
reduce the measurement of <5 by a factor of two, while a measurement of 623 will remain 
robust. In addition, it is possible to improve the existing upper limits on some of the NSI 
parameters by a factors of two. 
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FIG. 4. Effect of marginalizing over source and detector NSI parameters on precision measurements at 
ESSz^SB. Each panel shows the allowed region in the test 023 — ^ plane, when the NSI parameters are taken 
to be zero in the data. The red, green and blue curves represent the 68%, 90% and 95% C.L. contours, 
respectively. The solid contours show the effect of marginalization over the NSI parameters, whereas the 
dashed contours are for the standard oscillation scenario in the absence of NSIs. 
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FIG. 5. Effect of marginalizing over source and detector NSI parameters on precision measurements at 
ESSz^SB. Each panel shows the allowed region in the test 023 — S plane, when the NSI parameters are taken 
to be non-zero in the data. The true values of the amplitudes of the NSI parameters are assumed to be half 
of their 90% C.L. bounds from Ref. [2]|. The red, green and blue curves represent the 68%, 90% and 95% 
C.L. contours, respectively. The solid contours show the effect of marginalization over the NSI parameters, 
whereas the dashed contours are for the standard oscillation scenario in the absence of NSIs. 
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FIG. 6. Precision measurements at ESSj^SB for the case where NSIs are present in nature, but are not 
scanned for. Each panel shows the allowed region in the test 023 — S plane, when the NSI parameters are 
taken to be non-zero in the data. The true values of the amplitudes of the NSI parameters are assumed to 
be half of their 90% C.L. bounds from Ref. [2l|. The red, green and blue curves represent the 68%, 90% 
and 95% C.L. contours, respectively. The solid (dashed) contours show the allowed region without (with) 
marginalization over the NSI parameters. 
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FIG. 7. Limits on the amplitudes of the NSI parameters imposed by ESSi^SB data: as a function of the 

test value of the amplitudes of the NSI parameters, when the true value is zero. In each panel, all neutrino 
parameters (apart from the one indicated) have been marginalized over. The dotted lines from bottom to 
top show the 68%, 90% and 95% C.L., respectively. 
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FIG. 8. Precision on the amplitudes of the NSI parameters from ESS^SB data: as a function of the 

test value of the amplitudes of the NSI parameters, when the true value is non-zero. The true values of the 
amplitudes of the NSI parameters are assumed to be half of their 90% C.L. bounds from Ref. 2]|. In each 
panel, all neutrino parameters (apart from the one indicated) have been marginalized over. The dotted lines 
from bottom to top show the 68%, 90% and 95% C.L., respectively. 
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FIG. 9. Precision on the phases of the NSI parameters from ESSi/SB data: as a function of the test 

value of the NSI phases. The true values of the amplitudes of the NSI parameters are assumed to be half 
of their 90% C.L. bounds from Ref. 2l|, while the true values of the phases are taken to be zero. In each 
panel, all neutrino parameters (apart from the one indicated) have been marginalized over. The dotted lines 
from bottom to top show the 68%, 90% and 95% C.L., respectively. 
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FIG. 10. Role of the near detector in precision measurements at ESSz/SB. Each panel shows the allowed 


region in the test 023 ~ ^ plane, when the NSI parameters are taken to be non-zero in the data. The true 
values of the amplitudes of the NSI parameters are assumed to be half of their 90% C.L. bounds from 
Ref. [2l|. The red, green and blue curves represent the 68%, 90% and 95% C.L. contours, respectively. The 
solid (dashed) contours show the allowed region with (without) the near detector. 
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